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We numerically calculate, at the edge of chaos, the time evolution of the nonextensive 
, entropic form Sq = [1 — Pil/I^' ~ 1] (with Si = — ^^^pilnpi) for two families of 

one-dimensional dissipative maps, namely a logistic- and a periodic-like with arbitrary 
inflexion z at their maximum. At t = we choose initial conditions inside one of the 
W small windows in which the accessible phase space is partitioned; to neutralize large 
' fluctuations we conveniently average over a large amount of initial windows. We verify 

that one and only one value q* < 1 exists such that the limt^oo limvi^^oo limjv^oo Sq{t)/t 
is finite, thus generalizing the (ensemble version of) Kolmogorov-Sinai entropy (which 
corresponds to g* = 1 in the present formalism). This special, ^-dependent, value q* 
numerically coincides, for both families of maps and all z, with the one previously found 
through two other independent procedures (sensitivity to the initial conditions and mul- 



C/3 ' tifractal f{ct) function). 
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I. INTRODUCTION 

> 

. In the area of nonlinear dynamical systems, a number of studies have recently addressed the sensitivity 
to initial conditions, multifractality and the behavior of the Kolmogorov-Sinai (KS) entropy of these 
04 1 systems. In this context, it is worth mentioning that these attempts include both dissipative (low- 
8.nd high- [|| dimensional maps, symbolic sequences ^) and conservative (long-ranged many-body 
Hamiltonians Qjs], conservative maps [^) systems. Before describing the purpose of the present work, 
we need to briefly review the sensitivity to the initial conditions and the multifractality of the chaotic 
attractor. 

Firstly, let us focus on the concept of the sensitivity to initial conditions. As already well studied in the 
^ ^ literature for one-dimensional systems it is convenient to introduce ^(t) = lim^2,(g)_,Q -^f^ , where 

I ' Aa;(0) is the discrepancy of the initial conditions at time t — 0, and Ax{t) its time dependence. It can 
^ . be shown that ^ satisfies the differential equation d^/dt = Ai^, where Ai is the Lyapunov exponent, thus 
f (t) = exp(Aii). Consequently, if Ai < (Ai > 0) the system is said to be strongly insensitive (sensitive) 
to the initial conditions (and intermediate rounding). On the other hand, if Ai = 0, then the function ^ 
is expected to satisfy the diferential equation d^/dt — A^^', hence 



C(t) = [l + (l-g)A,t]^/(^-^' (ge7^), (1) 

5^ , which recovers the standard exponential case for q — 1, whereas 1 yields a power-law behaviour. 
If q > 1 (g < 1) and A^ < (A^ > 0) the system is said to be weakly insensitive (sensitive) to the 
initial conditions. Although asymptotic power-law sensitivity to the initial conditions has since long been 
observed ||lO|-p^, eq.(l) (which in fact corresponds to the power- law growth of the upper-hound of £,{€)) 
provides in principle a more complete description since t does not need to satisfy t » 1. At the onset 
of chaos (where it is Ai = 0), this upper bound (^ oc t^/^~'^) allows us to estimate the value q* of index q 
for the map under consideration. This method has been successfully used for a variety of maps such as 
logistic [0, z-logistic ||^, circle Q and z-circular [|| maps. 
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Secondly, we consider another interesting property of dynamical systems: the geometrical aspects of the 
attractor at the onset of chaos. In order to describe the scaling behavior of the critical dynamical attractor 
it is convenient to introduce a multifractal formalism . In this formalism, it is possible to introduce 

a partition function xq{-^) — ^iLiPf j where pi represents the probability on the ith box among the 
N boxes of the measure (it is necessary to warn the reader that we use here Q instead of the standard 
notation q of multifractal literature in order to avoid confusion with the present index q) . In the N ^ oo 
limit, the contribution to this partition function is proportional to N~^^^\ which comes from a subset 
of all possible boxes whose number scales as Nq oc N^'^'^\ where f{Q) is the fractal dimension of the 
subset. The content on each contributing box scales as Pq cx A^^"(<3) and all these exponents are related 
by a Legendre transformation t{Q) = Qa{Q) — f{Q)- The multifractal measure is then characterized by 
multifractal function /(a), which reflects the fractal dimension of the subset with singularity strength a. 
At the end points of the /(a) curve, this singularity strength is associated with the most concentrated 
[namely, amin = linig^oo ck{Q)] and the most rarefied [namely, amax = limg— »-oo o:{Q)] regions on the 
attractor. Recently, the scaling behaviour of these regions has been used to estimate the power-law 
divergence of nearby trajectories and a new scaling relation has been proposed [^ as 

111 

(2) 



1 Q ^min ^max 

This relation constitutes a completely different method for the calculation of the index q* . Previous 
works [|l|-Q have shown, for various map families, that the results of these two abovementioned methods 
of calculating the q* values are the same within a good precision. 

We are now prepared to describe the purpose of the present work, i.e., a specific generalization of the 
KS entropy Ki. For a chaotic dynamical system, one can define this entropy as the increase, per unit 
time, of the standard Boltzmann-Gibbs entropy 5*1 = — l^P*- Moreover, it is well-known that 

the KS entropy is deeply related to the Lyapunov exponents since the Pesin equality |l^ states that, for 
vast classes of nonlinear dynamical systems, Ki = Ai if Ai > and Ki = otherwise. Strictly speaking, 
as we shall detail later on, the KS entropy is defined in terms of a single trajectory in phase space, using a 
symbolic representation of the regions of a partitioned phase space. However, it appears that, in almost 
all cases, this definition can be equivalently replaced by one based on an ensemble of initial conditions. 
This is the version we use herein. 

The marginal cases, i.e., those for which Ai = 0, include period doubling and tangent bifurcations as 
well as the onset of chaos. For these cases, a generalized version of the KS entropy Kg has been introduced 
iQl as the increase rate of a proper nonextensive entropy 

= l^^EiMM! . (3) 

9-1 

This nonextensive entropy enables a generalization of Boltzmann-Gibbs statistical mechanics [|l6|; it 
recovers the standard Boltzmann-Gibbs entropy 5*1 — —"^^iPi^^Pi in the q — > 1 limit. A general 
review of the properties of this entropy and related subjects can be found in [|l^. So, for the generalized 
version of the KS entropy, it has been proposed [Q 

K„ = lim lim lim (4) 

^ t^ooW^oo N^oo t 

where t is the time steps, W is the number of regions in the partition of the phase space and N is the 
number of points that are evolving with time. The Pesin equality itself is expected to be generalizable 
as follows: Kq = Xq (or some appropriate average) if Ag > and Kq = otherwise. 

These ideas have been used very recently by Latora et al |l8j to construct a third method for the 
calculation of the q* values. They conjectured that (i) a special q* value exists such that Kq is finite 
for q — q* , vanishes for q > q* and diverges for q < q* , (ii) this value of q* coincides with that coming 
from other two methods of finding q* (namely, from eq.(l) and eq.(2)). Latora et al have checked these 
conjectures with numerical calculations for the standard logistic map and found that the growth of Sq{t) 
is linear when the value of index q equals q* ~ 0.2445, which strongly supports the point that all three 
methods yield one and the same special value of the index q* ! Although the results of Latora et al provide 
strong evidence in favor of this scenario, it is no doubt necessary to study more general maps along these 
lines in order to see whether the results of this method reproduce those of the previous two methods. 
The aim of this work is to address this question by studying the z-logistic and ^-periodic maps. 
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II. NUMERICAL RESULTS 



We first recall the z-logistic map 



Xt+i 1 - a\xt\ 



(5) 



where 1 < z, < a < 2, —1 < Xt < 1, and the z-periodic map 



Xt+l = d cos TT Xt - - 



z/2 



) 



(6) 



where 1 < z, < d < oo, —d < Xt < d. These maps display a period-doubling route to chaos and the 
parameter z is the inflexion of the map at its extremal point. Both maps belong to the same universality 
class, namely they share the same value q* for a given z. The values of flc and dc at the onset of chaos, as 
well as the q* values calculated from eq.(l) and eq.(2) are indicated in the Table for three representative 
values of z. 

Now we can describe the numerical procedure that we used as the third method for calculating q*. 
This method was first introduced in ^ for conservative systems and then used by Latora et al | [lq| for 
the standard {z = 2) logistic map. We partition the interval — 1 < a; < 1 {—d < x < d) into W equal 
windows for the z-logistic map (z-periodic map). Then we choose (randomly or not) one of these windows 
and select (randomly or uniformly) N initial values of x (all inside the chosen window) for the z-logistic 
(z-periodic) map at a given value of a {d). As t evolves, these N points typically spread within the 
interval of phase space and this gives us a set {Ni{t)} with — which consequently 

yields a set of probabilities {pi{t) = Ni{t)/N}. Since for i = all TV points are inside the chosen window, 
5,(0) = and, as time evolves, Sq{t) starts increasing, although bounded by the equiprobability value 
(i.e., by {W-^'"^ — — q) for q I, and InW for q = 1). Let us anticipate our main result. For N 

and W increasingly large (and always satisfying N » W) we observe, in all cases, that a /mear increase 
of Sq with time emerges only for a special value of q, namely q*{z). We define the corresponding slope 
as the generalized Kolmogorov-Sinai entropy Kg. For values of a and d for which Ai > we verify that 
q* — I. But at the edge of chaos, we obtain q* < 1. 

In contrast with the cases where strong chaos exists, considerably large fluctuations appear in Sq{t) at 
the chaos threshold, due to the fact that the attractor occupies only a tiny part of the available phase 
space. These fluctuations make uneasy the task of determining the value of q for which linearity is present. 
In order to neutralize these fluctuations, we adopt a procedure of averaging over the "efficient" initial 
conditions as introduced in [|l8|. The procedure is as follows: we choose the initial distribution of N 
points in one of the W cells of the partition and count the total number of occupied cells as time evolves 
between say t = I and t = 50; this number is a measure of the efficiency of that particular window in 
spreading itself. After studying each one of the W/2 cells in the interval < a; < 1 (0 < cc < d) for the 
z-logistic (z-periodic) maps, the averaging is done over the cells for which the total number of occupied 
cells is larger than a fixed threshold, say 5000 (see Fig. 1). 

Our results for typical values of z for the z-logistic and the z-periodic maps at the threshold of chaos 
are indicated in Figs. 2-4. We used iV = 10 x = 10^ for the z-logistic maps, and A^=10xW^ = 6x 10^ 
for the z-periodic maps. The values we have used are large enough so that the time evolution of the 
entropy in the intermediate time region (after the initial transient and before the saturation) does not 
depend much on W. This point has been illustrated in Fig. 2, where we plot the time evolution of Sq{t) 
for two values of W for the case of the z = 1.75 logistic map. It is clear from this figure that the starting 
points of saturation of the curves shift to larger time as W increases, but in the intermediate time region 
the curves are almost insensitive to W. 

In all cases, the growth of Sq{t) in this intermediate time region is found to be linear when q — q* , 
whereas it curves upwards for q < q* and downwards for q > q*, i.e., Kq is finite for q = q*, diverges 
for q < q* and vanishes for q > q*. In order to provide quantitative support to this behavior, we fit the 
curves with the polynomial S{t) — A + Bt + Ct^ in the interval [ti,i2] characterizing the intermediate 
region. We define the nonlinearity coefficient R = C{ti + t2)/B as a measure of the importance of the 
nonlinear term in the fit; R vanishes for a strictly linear fit (see the Insets of Figs. 2-4). The times ti 
and t2 are respectively defined as the end of the initial transient (during which 5*^ is roughly constant) 
and the beginning of the saturation; R{q) is almost independent from (ti, for any set of specific values 
(^1,^2) which satisfy the above characterizations. We also see in the insets that the signs of R on both 
sides of q* are consistent with the curvatures of Sq{t) in the intermediate region. 
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III. CONCLUSIONS 



Let us make a few general considerations in order to better expose the meaning of the present effort 
and its results. The Kolmogorov-Sinai entropy Ki is an important concept in chaotic dynamical systems, 
either dissipative (like the one-dimensional maps considered here) or conservative (like classical many- 
body Hamiltonians, satisfying the Liouville theorem) . Its definition is based on a partition of the accessible 
phase space in a set of W subspaces that are visited, along time, in some complex order starting from a 
single initial point in phase space. If we associate to each subspace a symbol, we shall have possible 
words of length Af. These W-^ words are visited along time with probabilities {tt;} {I = 1,2, ...,W^). 
This set of probabilities enable the calculation of the Boltzmann-Gibbs-Shannon entropy Si{{tti}) — 
— ^jTTilnTT;. In the limit Af — > oo, S'i({7r;}) is proportional to J\f if the system mixes exponentially 
quickly (i.e., positive Lyapunov exponents), and Ki is defined as the supremum of the lim7v^_+oo Si{Af)/Af. 
Although we are not aware of any general proof, a common belief exists that this computationally quite 
heavy definition can be conveniently replaced by the one we have used here, based on an ensemble of 
initial conditions instead of on single trajectories. This is to say, we choose N initial conditions within 
one of the W subspaces, and we follow along time the set of probabilities {pi} {i = 1, 2, W) associated 
with the occupancies of those subspaces. The set {pi} enables the calculation of Si{{pi}) — — pi Inp^, 
and K\ is expected to be the supremum of the limt^oo •S'l (i) jt ■ 

The scenario we have just described is the standard one. What we have exhibited here is a step forward, 
although only for very simple cases, namely the maps herein considered. By analyzing the generalization 
[Ql Kq of the KS entropy based on the nonextensive entropy 5'g, either in its trajectory formulation or 
in its ensemble one, we exhibit a path that might be applicable to situations much more general than 
the maps focused on here. We have (numerically) shown that an unique value q* exists such that Kq 
is finite, being zero for q > q* and infinite for q < q*. If the nonlinear dynamical system is strongly 
chaotic (exponential mixing in phase space), then q* — 1, thus recovering the usual scheme. But if 
the system is only weakly chaotic (power-law mixing, i.e., vanishing Lyapunov exponents, but positive 
generalized Lyapunov exponents, like A^), then one essentially expects q* < 1. Let us emphasize that 
in the present work we have only addressed q* and not Kq* (nor A,*). Indeed, the averaging procedure 
we have implemented is expected to preserve linearity with time, but certainly not the supremum, or 
whatever analogous to that. Further analysis is certainly welcome. 

As in the case discussed by Latora et al we have verified that the value q* that has emerged is 
precisely the same previously obtained through two completely different procedures, namely the power- 
law sensitivity to the initial conditions (see eq. 1), and the multifractal structure of the chaotic attractor 
(see eq. 2). This uniqueness of the special value q* clearly provides to the whole scenario a very robust 
consistency. 

Let us end with a speculative question quite analogous to the one we have addressed here. If we 
consider two systems A and B that are independent in the sense that pfj^^ = pfpf, then we straightfor- 
wardly verify that Sq(A + B) = Sq(A) + Sq(B) + (1 — q)Sq(A)Sq(B). In other words, we have extensivity, 
superextensivity and subextensivity for q = 1, q < 1 and q > 1 respectively. Let us now assume that 
we have a classical d-dimensional many-body Hamiltonian system whose particles interact through two- 
body interactions which are nowhere singular and which (attractively) decay like at long distances. 
It is known that such a system is thermodynamically extensive ii a/d > 1 (short-range interaction) and 
nonextensive if < a/d < 1 (long-range interaction) |l7[ . For a/d > 1, if the system is thermodynam- 
ically large (say it contains M particles with M ~ 10^, any two, also thermodynamically large, of its 
subsystems can be considered as probabilistically independent in the above sense. Then, because of its 
well known extensivity, the special value of q to be associated clearly is q = 1. Consequently the entropy 
is generically expected to yield a finite value for limf^oo hniA/^oo Si(t, M)/M. The question we wish 
to leave open at the present stage is: for < a/d < 1, is there a special (possibly unique) value of q such 
that limt^oo limM-»oo Sq(t, M)/M also is finite? 
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Figure and Table Captions 

Figure 1 - Total number of occupied cells as function of the window rank (varying from rank 30000 to 
rank 60000, respectively corresponding to a; = and x = 1) for the z = 3 periodic map with W = 6x 10^. 
The horizontal line indicates the cutoff at 5000. Average for Sg{t) is done over the windows associated 
with values above the cutoff. 

Figure 2 - Time evolution of Sq{t) for two different values of W for the z = 1.75 logistic map. Inset: 
The nonlinearity coefficient R versus q (see text for details). 

Figure 3 - Time evolution of Sq{t) for three different values of q and W = 10^ for the z = 3 logistic 
map. Inset: The nonlinearity coefficient R versus q (see text for details). 

Figure 4 - Time evolution of Sq{t) for three different values of q and W = 6 x lO** for the z-periodic 
maps: (a) z = 1.75; (b) z = 2; (c) z = 3. Insets: The nonlinearity coefficient R versus q (see text for 
details) . 

Table - The values of Uc, dc and q* at the chaos threshold for typical values of the inffexion parameter 

z. 



Table 





z-logistic map 


z-periodic map 




z 






q* 


1.75 


1.35506075... 


0.77946132... 


0.100 


2.0 


1.40115518... 


0.86557926... 


0.2445 


3.0 


1.52187879... 


1.07848805... 


0.472 
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